Minimal seesaw model with tri/bi-maximal mixing and leptogenesis 
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We examine minimal seesaw mechanism in which the masses of light neutrinos are described with 
tri/bi-maximal mixing in the basis where the charged- lepton Yukawa matrix and heavy Majorana 
neutrino mass matrix are diagonal. We search for all possible Dirac mass textures which contain at 
least one zero entry in 3 x 2 matrix and evaluate the corresponding lepton asymmetries. We present 
the baryon asymmetry in terms of a single low energy unknown, a Majorana CP phase to be clued 
from neutrinoless double beta decay. 
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I. INTRODUCTION 

^sO ' Thanks to the accumulating data from experiments on the atmospheric and solar neutrinos experiments 0, 0> @] j 
we are now convinced that neutrinos oscillate among three active neutrinos. Interpreting each experiment in terms 
of two-flavor mixing, the mixing angle for the oscillation of atmospheric neutrinos is understood to be maximal or 
nearly maximal: sin 2 28 a t m — 1, , whereas the one for the oscillation of solar neutrinos is not maximal but large: 
sin 2 9 so i ~ 0.3, 0. The upper bound for Qrgaa sin# reQC < 0.2, was obtained from the non-observation of the 
disappearance of in the Chooz experiment Q with Am 2 < 10~ 3 eV 2 . The unitary mixing matrix is defined via 
v a = U a jVj (a = e, [i, r), where v a is a flavor eigenstate and t/j is a mass eigenstate. Including data from SNO0 

and KamLand 5] , the range of the magnitude of the MNS mixing matrix is given by 0, 0- IH] , 

/ 0.79 - 0.86 0.50- 0.61 
\U\= 0.24- 0.52 0.44- 0.69 0.63- 0.79 | (1) 
' \ 0.26 - 0.52 0.47- 0.71 

O ,' at the 90% confidence level. The existing data also show that the neutrino mass squared differences which induce 
the solar and atmospheric neutrino oscillations are Am 2 o/ ~ (7^2°) x ^ b eV 2 and Am 2 tm ~ (2.5tJ;g) x 10 -3 eV 2 , 
1 respectively. It can be readily recognized that the central values of elements in the mixing matrix in Eq.Q are 
^ , pointing an elegant form, which is called tri/bi-maximal mixing (llj. 

X : ( J* 

Utb = 73 ~72 I ' ( 2 ) 

\ L J_ i 

which consists of sm9 so i — -4= and sin# atm = -4=. There are some literatures [l^ which proposed textures of the mass 
matrix based on the particular mixing type Utb- 

On the other hand, the baryon density of our universe Vlsh 2 = 0.0224 ± 0.0009 implied by WMAP (Wilkinson 
Microwave Anisotropy Probe) data indicates the observed baryon asymmetry in the Universe |13t . 

V C b MB = = (6.5±°;4) x 10- 10 , (3) 

where ub , rig and n 7 are number density of baryon, anti-baryon and photon, respectively. The leptogenesis |l5j has 
become a compelling theory to explain the observed baryon asymmetry in the universe, due to increasing reliance on 
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the seesaw mechanism from experiments. Theory for lepton asymmetry requires two heavy right-handed neutrinos or 
more. For that reason, a class of models with two heavy right-handed neutrinos and 3x2 neutrino Dirac mass matrix is 
called the minimal neutrino seesaw models(MNSMs) which were intensively studied by several authors recently pjjll7T |. 
especially for simple Dirac mass textures that make prediction compatible with solar and atmospheric neutrino data. 

The main framework of our work is seesaw mechanism in bottom-up approach. We launch our analysis by taking 
Utb for mixing of light neutrinos and then investigate the structure of 3 x 2 Dirac matrix. That is, our concern remains 
on the combination of tri/bi-maximal mixing and MNSMs, in order to study the phenomenological implication of the 
high energy theory based on the low energy theory. One advantage of our framework is that physical observables 
can be explained in minimal terms of physical parameters. In Section II, we present the light neutrino mass matrix 
in terms of the mixing given in Eq.@ and mass squared differences measured in experiment. The mass matrix 
reconstructed in that way will constrain Dirac mass matrix. In subsections, depending on the type of mass hierarchy, 
possible 3x2 Dirac matrices will be examined carefully. In Section III, leptogenesis will be discussed in details based 
on the Dirac matrices investigated before. In section IV, we will present numerical results on leptogenesis in our 
scheme and a relationship between leptogenesis and neutrinoless double beta decay as well as the lower bound of Mi 
will be discussed focusing on the contribution from a single Majorana phase. 



II. DIRAC MASS MATRICES IN MINIMAL SEESAW 



In general, a unitary mixing matrix for 3 generations of neutrinos is given by 

U = R(e 23 )R(6 13 ,8)R(e 12 )P(<p,ip') (4) 

where i?'s are rotations with three angles and a Dirac phase 6 and the P = Diag ^1, e lv ^ 2 , e %v ^ 2 ^j with Ma- 
jorana phases ip and iff is a diagonal phase transformation. The mass matrix of light neutrinos is given by 
M v = UDiag(mi, rri2, m,3)U T , where mi, m 2 , m 3 are real positive masses of light neutrinos. Or the Majorana phases 
can be embedded in the diagonal mass matrix such that 

M v = UDiag(mi,fh,2,fh3)U T , (5) 

where U = UP^ 1 and m 2 = m 2 e ¥ and m 3 = m 3 e l,p . 

If the Utb in Eq.© is adopted for the U in Eq.JSJ), the light neutrino mass is 

1 

M v = m 1 | 1 | + | 1 1 1 | + ■^±_z±. | o 1 -1 | , (6) 

1 

which orients toward a minimal model of neutrino sector by removing an angle and the Dirac phase. With SNO 
and KamLand, data have narrowed down the possible mass spectrum of neutrinos into two types, Normal Hierarchy 
(NH), mi < m 2 < m 3 , and Inverse Hierarchy (IH), m 3 < mi < m 2 for MSW LMA. Those two types include the 
possibility of zero mass for a neutrino, which is necessarily followed by relegating one of the Majorana phases to the 
unphysical. In other words, the minimal model with the physical observables which the present experimental data 
guarantee can be obtained by Utb and dictating zero mass to one generation of neutrinos, where the non-zero physical 
parameters in the model consist of 2 masses, 2 angles, one Majorana phase. 

When only two of three neutrinos are massive, by accommodating the experimental results Am 2 ol = m 2 — m 2 
and Am 2 tm = |mj — m 2 | to the two types of mass hierarchies, one can obtain the following expressions for mass 
eigenvalues, 

mi = 




m 2 = y /Am 2 ol for NH (7) 

m 3 = v 7 Am L + Am atm 

m i = \/ ^mj tm - Amj ol 

m 2 = VAmL; for IH. (8) 

m 3 = 

Phase transformation P = Diag (l, e up / 2 , l) now can replace the phase transformation in Eq.(0} without loss of 
generality, whether NH or IH, so that one can single m 2 out in order to investigate the CP violating contribution of 
the Majorana phase. 
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Effective neutrino mass models with one zero mass eigenvalue involved in three active neutrinos can be generated 
naturally from MNSMs. In the basis the mass matrix Mr of right-handed neutrinos N R — (Ni,N 2 ) is diagonal, the 
model is given 



C = -I l M l Ir - T7Zm D N R + -N C R M R N R + h.c. 



(9) 



from which the light masses are derived through the seesaw mechanism, M„ = —mjjM^rrVjj in top-down approach. 
On the other hand, the matrix m_o is found as the solution to the seesaw mechanism in bottom- up approach once one 
launches the analysis with the light neutrino masses M„. Let M\ and M2 be the masses of right-handed neutrinos 
and Mij the elements of the matrix M v . The Dirac matrix, 



m D 



/M~[a 2 y/M 2 ~b 2 
/Mia 3 



(10) 



/NUb 



2 "3 



is resulted in with 



M n -b(, h = JMx 



1 



1 



aiM H - (TihJMnM, 



M U M 23 = 



61 M u + amJMiiMa - Mfi 



(11) 



M12M13 + <7 2 <7 3 J (M11M22 - Mf 2 ) (M11M33 - Mfg) 



where the i is 2 or 3, the <7j is a sign ±, and the sign of a% is fixed as positive. The solution in eq. (|12|l was first 
derived and formulated in Ref.^J]. It is clear that only 5 parameters out of 6; a\, 61, ai, bi, can be specified in terms 
of the elements of M„. There are various ways to decrease the number of parameters in Dirac matrix, posing one or 
more zeros or posing equalities between elements for the matrix texture. It is known that texture zeros or equalities 
among matrix entries can be generated by imposing additional symmetries to the theory. 

In this paper, we focus on posing one or more zeros in Dirac matrix, and show that only one-zero textures are 
allowed for NH and only one-zero and two-zero textures are allowed for IH, accompanied with the low energy mixing 
Utb- On the other hand, from Ea. l|llf> . one can recognize that, if there exists a kind of symmetry between entries 
in M„, the Dirac matrix also has a symmetry in certain entries inherited from the symmetry of the M v . So, there 
are a number of patterns with equalities among the entries in Dirac matrices obtained with one or two zeros, as a 
consequence of maximal mixing. 



A. Normal Hierarchy 

With mi = 0, the neutrino mass M v in Eq.© reduces to 

M„ = 




(12) 



where d — 3m3/2m2, which, using Eq. l|ll|) . gives rise to Dirac matrix with the following entries: 



ai = \/m 2 /3 - b\, bi = \Jm 2 lZ - a\, 

a>i = o>\ — Cj&iVa, bi = b\ + atai^/d, i = 2, 3 



(13) 



where a 2 oz = —1. Depending on the position of texture zero, the types of Dirac matrices can be classified as follows; 



• NH 1-a : 61 = 0, ai = a i = a 2 = a 3 , b 2 = -63 = a 2 ^m 3 /2 

• NH 1-b : cti = 0, ax <-> b%, ai <-> bi in NH 1-a 

• NH 2-a : b 2 = 0, a 2 = ^rh 2 /3 + m 3 /2, ai = o 2 /(l + d), a 3 = a 2 (l - d)/(l + d), 61/ai = -a 2 Vd, b 3 = 26i 
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• NH 2-b : a 2 = 0, ax <-> bt, a, L <-> bi in TVii 2-a 

The matrix in Ea. (|12|l features the equalities between M22 and M33 and between M\i and Mi 3 as a consequence of 
the maximal mixing of atmospheric neutrinos. In the case with 62 = 0, it can be recognized that the ratios of a\ to 0,2 
and a-i to 02 inherit those of Mi 2 to M 2 2 and M23 to M 2 2, respectively. 63 = or 03 = cases will not be presented 
as an independent case, since it can be made by exchanging 63 with 62 and 03 with 02 from NH 2-a and 2-b. 



B. Inverse Hierarchy 

With m 3 = 0, 



2 x- 1 03-1 
x + 1/ 
1+1/2 

where x = 1712/1121, which, using Ea. Hlljl . gives rise to Dirac matrix with the following entries 



Mu = IT ' •'' + 1/2 ■'' + :l/2 I 1141 



ai = y/m\{x + 2)/3 — 6f , 61 = y/m\{x + 2)/3 — of, 

Oi = ((x- l)ai -3o-6i i/xT^ j /(« + 2), h = ((x - l)h + 3aa ly /x/2) j(x + 2), i = 2, 3 



(15) 



where cr = cr,; and (T2O3 = 1. The equality M22 = M23 = M33 in Ea. (|14() . which is again a consequence of the 
maximal mixing of atmospheric neutrinos, constrains the elements of the Dirac matrix such that a 2 = a 3 , £> 2 — &3- 
Hence, texture with a single zero included appears only if a\ = or b\ = 0, while texture with two zeros appears if 
0-2 — a-3 — or 62 = &3 = . 



• IH 1-a : 61 = 0, a% = \/rri2/3 + 2mi/3, a, = ai(x — l)/(a; + 2), fo/a.; = 3<jy/x/2/(x — 1) 

• /if 1-b : a\ = 0, ai <-> &i, «-> 6^ in Jii 1-a 



• JiJ 2-a : 62 = &3 = 0, ei2 = 03 = \Jm2l3 + mi/6, ai = 2a2(x — l)/(2x + 1), 61/ai = —3ay/x/2/ (x — 1) 

• Jii 2-b : <22 = 03 = 0, ai «-> 61, a% «-> &i in JiJ 2-a 

Listed are all the cases with one or more texture zeros in Dirac matrix derivable from the light neutrino mass with 
Utb, whether NH or IH. In the following, the eligibility of each case to generate the CP asymmetry for leptogenesis 
will be examined. 



III. LEPTOGENESIS 



The baryon asymmetry Eq.© can be rephrased 

Y B = UB ~ n * ~ (8.8 - 9.8) X HP 11 . (16) 
s 

The n 7 is the photon number density and the s is entropy density so that the number density with respect to a co- 
moving volume element is taken into account. The baryon asymmetry produced through sphaleron process is related 
to the lepton asymmetry 0, ^| by 

Y B = oYb—l = — -t^l, (17) 
a — 1 

where a = {%N F + AN H )/{22N F + IZN H ), for example, a = 28/79 for the Standard Model(SM) with three generations 
of fermions and a single Higgs doublet, N F = 3, Nh = L The purpose of this work is to estimate whether the Yukawa 
interaction which produces the light neutrinos with the mixing Eq. through the seesaw mechanism can also generate 
a sufficient lepton asymmetry for the observed baryon asymmetry. The generation of a lepton asymmetry requires 
the CP-asymmetry and out-of-equilibrium condition. The Yl is explicitly parameterized by two factors, e, the size of 
CP asymmetry, and k, the dilution factor from washout process. 



(18) 
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where g* ~ 110 is the number of relativistic degree of freedom. The e, is the magnitude of CP asymmetry in decays 
of heavy Majorana neutrinos pcLISl) . 

TjN, -> IB) - T(N t -> IB*) 
£i r(iVj — » + T(Ni — > ' 1 j 

where i denotes a generation. When one of two generations of right neutrinos has a mass far below that for the other 
generation, i.e., Mi < M 2 , the in Eq.JT^Jl is obtained from the decay of Mi p2 l |23L |24| . 

1 Im [(mjtmj)^] /MA 
El 8™ 2 (m D t TOD ) n 1 \ Ml ) ' lZUJ 

where v = 174 GeV and / (M2/M1) represents loop contribution to the decay width from vertex and self energy and 
is given by 



f(x) = x 



1 - (1 + .x 2 ) In — + 



(21) 



for the Standard Model. For large value of x, the leading order of f(x) is (— 3/2)x 1 . 

It is convenient to consider separately the factor that depends on Dirac matrix in E\ in Ea. (|20[) at this stage. 



Im [{m D ^m D )\ 2 ] Im [(a x bi + a 2 b 2 + a 3 6 3 ) _ A / 00 n 

— 7 1 \ = M 2 i — I2~ri — 12TI — 12 = M 2 Ai, (22) 

(to d Tto d )h | ai |^ + | ffl2 |-J + | a3 |^ 

where a's and 6's arc defined in Eq.(ll). From a number of types of matrices with a texture zero derived in Ea. (jl^f) 
and Ea. (|15J) . only 6 different non-zero values of Ai's can be evaluated. Those particular Dirac matrices to contribute 
the imaginary parts are the matrix with b 2 = and that with a 2 = for NB, and the matrix with b\ = 0, that with 
a\ = 0, that with b 2 = 63 = 0, and that with a 2 = 03 = for IB. For NB, if ai = 0, or 61 = 0, the (m D m,D)i2 
vanishes from the trivial relation between entries. Applying Ea. i|13[l and Eq. l|15[) for Ea. H22|) . one can find that each 
type of Dirac matrix gives rise to Ai as follows; 

A 1 (NB 2 - a) = 6to,2TO 3 (to 2 ~ m 2 ) sin tp (23) 

(2m,2 + 3m§)v 4»^2 + ® m 3 + 12m 2 m 3 cos tp 

« ,„ TTT „ ,\ — 6(mo — to 2 ,) sin tp , 

A 1 (NB2-b) = — v 3 21 = , (24) 

5y/ 4m 2 , + 9to 3 + 12rn2TO 3 cos tp 

where m 2 and m 3 are given in terms of Ato 2 o; and Am 2 tm in Eq.JJJ, 

A 1 (IB 1 - a) = 2ui\ m 2 (to 2 to 2 ) sin tp ^ 

(2m\ + to 2 .) \J 4m 2 + TO i + 4toito 2 cos tp 

A 1 (IB 1 - b) = 2{ml-ml)s^tp (2g) 

3y 4m 2 + to 2 + Am\m 2 cos tp 

-2toito 2 (to| -to 2 ) sin y 
Ai(IB2-a) = — — - (27) 

(TO^ + 2TO-2 ) y TO I + 4TO-2 + 4TO1 TO2 cos tp 

A l{ IB 2 - b) = -pjM^4^ , (28) 
3y toj + 4to2 + 4toiTO2 cos </? 

where toi and TO2 are given in terms of Am 2 sol and ATO 2 tm in Eq.JSJ. Thus, for M 2 3> Mi case, the CP asymmetry 
in Ea. l|20|) reduces to e\ « 3/(167rw 2 )Af 1 A 1 , which is now parameterized by the lightest mass of heavy neutrino M\ 
and Majorana phase tp. The sign of e\ depends on the position of a texture zero in a row of Dirac matrix. 

The k in Ea. (|18f) is determined by solving the full Boltzmann equations. The k can be simply parameterized in 
terms of K defined as the ratio of Fx the tree-level decay width of N\ to B the Hubble parameter at temperature 
Mi, where K = T\/B < 1 describes processes out of thermal equilibrium and k < 1 describes washout effect |Tflj . 

<U for 10 < K < 10 6 , (29) 



if (In if) ' 6 



- 1 for < K < 10. (30) 
2V^ 2 + 9 V ' 



6 



The decay width of Ni by the Yukawa interaction at tree level and Hubble parameter in terms of temperature T 
and the Planck scale M p i are Y\ — (rn^rriDJ Mi/ (8irv 2 ) and H = 1.66gl^ 2 T 2 /M p i, respectively. At temperature 



T = Mi, the ratio K is 



K _ m pi (w^mg)ii 

1.66^(8^2) Mi ' 1 ' 

which reduces to, using the Dirac matrices in Eq. l|lU|) . 

K *J^y (|ai| 2 + |a 2 | 2 + |a 3 | 2 ), (32) 

where all fixed numbers are included in a factor of order. As done for the Ai's, one can apply Eq.(JT3} and Ea. l|15|) 
for Eq. (l: j2l to find dilution factor re when the decay width is determined by Yukawa couplings in each type of Dirac 
matrix. For the six types of Dirac matrices that are eligible for the CP asymmetry as in Eas. (|23|) - (|28|) . the ratio K 
for each case is 

( 2m?, + 3m?, 5m,2m,3 ) , . 

K(NH 2 -a, -b) w i 2 3 ; (33) 

(10 _3 eV) i/4m2 + 9m| + 12m2m3 cost/3 

( 2m? + m?, 3mims> ) 

K(IH 1 - a, -b) « i 1 2 ' ' (34) 

(10 3 eV)\/4m? 1 + m^ + 4niini2 cos<^ 

( m? + 2m?, 3mim9 ) 

K{IH2-a, -b) « 11 2 ' ' =, (35) 

(10 3 eV)-\/m? 1 + 4m| + 4mim2 costp 

which shows that the dilution factor also depends on the phase <p, but it does not significantly affect the order of 
magnitude. Out of all the types of Dirac matrices examined, there is no such a case that Yukawa couplings originate 
decays of neutrinos N\ which satisfy the out-of-equilibrium condition K < 1 at T = Mi. The washout effect of 
asymmetry is most suppressed with the Dirac matrix of type NH 2 — b, where, depending on ip, the dilution factor 
ranges from 0.010 to 0.013, the amount of asymmetry survived from washout is at most about 1%. When T < Mi, 
the Boltzmann equations still depict the finite value of re as M%/T increases for the universe evolution 



IV. DISCUSSION 



Based on the formulation of the leptogenesis derived in the previous section, we numerically analyze baryon asym- 
metry for each case classified as NH or IH. For the numerical calculation, we take Am 2 oi = 7.0 x 10~ 5 eV 2 and 
Am 2 4m = 2.5 x 10~ 3 eV 2 as inputs. 

Consider a model with neutrino masses in normal hierarchy. In Fig.l, we plot the baryon asymmetry Yg as a 
function of the Majorana phase ip for NH2 — b. The different curves correspond to Mi = 2.0 x 10 11 to 2.0 x 10 13 
GeV for fixed M2/M1 = 5. We note that we can choose any reasonable M2/M1 value which can protect L- violating 
processes with Ni from the wash-out when T < M2 ■ As expected from Eq. I|20|) , the value of Yb for a fixed <p increases 
with Mi. The horizontal line in Fig.l presents the current cosmological observation of Yb given in Ea. (|16(l . From 
the analysis, we see that the current observation on Yb constrains the lower bound of Mi, which turns out to be 
Mi > 2.0 x 10 11 GeV. It is clear that the CP asymmetry in high energy is almost proportional to the imaginary 
part of Majorana CP contribution in low energy from Eas.l|22 |) -(|28 |) . Thus, the plots show that the lower bound of 
Mi to generate the observed baryon asymmetry should be raised if the imaginary contribution of low energy phase 
is decreased as the ip approaches or it. In all aspects of the prediction of Yb, NH2 — a and NH2 — b are quite 
similar to each other except an overall factor. The Yb for NH2 — b is enhanced from both the enhancement of CP 
asymmetry, Ai(6)/Ai(a) ~ 3.6, and the suppression of wash-out effect, re(6)/re(a) ~ 4.5. The lower bound of Mi with 
<p = tt/2 is pulled down to 2.0 x 10 11 GeV for NH2 - b, whereas that for NH2 - a is 3.2 x 10 12 GeV. Suppressing a 
certain Yukawa coupling by putting a texture zero can vary the amount of the asymmetry by order of magnitude. 

In Fig. 2, we plot Yb as a function of the Majorana phase <p for I HI — a. The different curves correspond to 
Mi = 5.5 x 10 13 to 5.5 x 10 15 GeV for fixed M 2 /Ml = 5. As in NH, we obtain a lower bound on Mi > 5.5 x 10 12 
GeV for IH. The prediction of Yb for IH with the same value of Mi is smaller than that for NH because Ai for IH 
is proportional to m 2 , — m 2 which corresponds to the solar mass squared difference, while Ai for NH is proportional 
to m 2 ; — mf, which corresponds to the atmospheric mass squared difference. We expect from Eqs. (25-28) that the 
predictions of Yb's for other cases of IH are almost the same as that for IH 1-a because mim,2/(2m 2 + m 2 ,) ~ 1/3. 
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Although the Majorana phase is not detectable through neutrino oscillations, it may affect the amplitude of neu- 
trinolcss double beta decay. Thus, one can anticipate that there may exist a correlation between leptogenesis and 
neutrinoless double beta decay in our scenario where the heavy Majorana neutrino mass matrix and the charged-lepton 
Yukawa matrix are both diagonal. The neutrinoless double beta decay amplitude is proportional to the effective Ma- 
jorana mass |(m ee )| which can be written in the form: 

l(mee)| = |5»ie^| 

i=l 

= I ^' ( NH) , m) 

\ I (4m? + m| + 4toito 2 cos^) 1 ' 2 , (IH) v ; 

where (pi are Majorana CP-violating phases. The |(m ee )| depends on the CP phase (p only with inverted hierarchy, so 
that one can draw a simple correlation between leptogenesis and neutrinoless double beta decay only for the particular 
case. In Fig. 3, we present a correlation between Yb and |(m ee )| for IH 1-a. The inputs are taken to be the same 
as in Fig. 2. As the value of |(m ee )| approaches to that with <p = ir/2, the asymmetry is enhanced and the bound of 
Mi becomes lower. The lower bound of Mi as a function of Majorana phase or that of effective Majorana mass is 
obtained from the current cosmological observation of Yb- In Fig. 4, we present a correlation between the lower bound 
of Mi and |(m ee )|. 

We examined the minimal seesaw mechanism of 3 x 2 Dirac matrix by starting our analysis with the masses of 
light neutrinos with tri/bi- maximal mixing in the basis where the charged-lepton Yukawa matrix and heavy Majorana 
neutrino mass matrix are diagonal. We found all possible Dirac mass textures which contain one zero entry or two in 
the matrix and evaluated the corresponding lepton asymmetries. The baryon asymmetry can be presented in terms 
of low energy observables, where only one Majorana CP phase among them remains yet unknown. The numerical 
work exhibits the dependence of both the size of baryon asymmetry and the lower bound of Mi upon the low energy 
CP phase to be clued from neutrinoless double beta decay. 

Note added : After completing this wrok, we have been noticed that similar analysis for the hierarchical case in 
supersymmetric seesaw model appeared in Ref. |27j . 
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FIG. 1: Yb as a function of Majorana CP phase for case NH 2-b, with various values of Mi where M2/M1 = 5. The horizontal 
lines are the current cosmological bound of Yb ■ 
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FIG. 2: Yb as a function of Majorana CP phase for case IH 1-a, with various values of Mi where M2/M1 — 5 



1e-06 
1e-07 
1e-08 |r 
1e-09 
1e-10 
1e-11 
1e-12 r 
1e-13 



n rrz r 

5.5x 10 ;GeV - 
5.5 x 10 *GeV - 
5.5 x 10 15 GeV 




0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 

<m ee > eV 



FIG. 3: Y B as a function of |(m ee )| for case IH 1-a, with various values of Mi where M 2 /Mi = 5 
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FIG. 4: The lower bound of Mi as a function of |{m ee )| for case IH 1-a 



